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Abstract. Semi-analytic models of X-ray clusters and groups of galaxies, relying on the idea that there was a non- 
gravitational energy injection in these systems, are able to reproduce many observed correlations, in particular 
the Lx — T relation and the "central entropy floor" in groups. Limiting models exist which describe the behaviour 
of clusters and groups separately, but no analytic modeling has yet been found to unify both mass ranges. It is 
the aim of this paper to provide such an analytic model. Our descri ption relies on a now standard description of 
the shock thought to occur in these systems near the virial radius (Cavaliere et al., 1998), the isothermality and 
spherical symmetry of the intracluster medium, as well as the reinterpretation of observed quantities (like the 
X-ray luminosity, the gas mass Mjcm and the central SZ effect yo) in terms of the specific entropy. This allows 
the derivation of analytic expressions for several observed correlations [Lx — T, Mjcm —T,yo — T,...) and their 
normalisation encompassing both the group and the cluster regimes. The analytic predictions compare very well 
with observations, as well as with more elaborated semi-analytic schemes. This agreement allows a reinterpretation 
of the Lx — T relation (via the quantity Lx/T"^^^) and the yo —T relation (via yo/T^^^) as indirect measures of 
the non-gravitational entropy content of groups and clusters of galaxies. Finally, we emphasize the need for shock 
heating, even in the group mass range : shocks can not be completely supressed in groups (and thus groups can 
not be entirely isentropic) unless an unacceptably high entropy floor is needed in order to break the self-similarity 
in the Lx — T relation. 



Key words. Hydrodynamics 
: analytical 



Shock waves - clusters : galaxies : general - X-rays : galaxies : clusters - Methods 



1. Introduction 

For ten years, it has been known that X-ray observations 
of clusters of galaxies can not be reproduced in simple self- 
similar models, where the central gas density is propor- 



tional to dark matter density. Kaiser (1991 ) and Evrard 



& Hen 



active galactic nuclei. Several semi-analytic schemes that 
elaborate on this idea have been proposed to reproduce 
the curvature of the Lx — T relation in the groups' mass 
range (where the eff ects of the preheating a re thought to 
be the highest, see Cavaliere et al., 1997; Bower, 1997 



y (1991 ) advocaLed a piehcaLlug uf LliLi gas befui 



Valageas fc Silk, 1999; Valageas fc Schaeffer, 200C; Wu 



et al., 2000| ; p3ower et al., 2000 ). In particular, [Cavaliere 



it fell 



iiLu Lhc clusLei 



Llic ubseived 



lusLei puLeuLial Lu lecuvei 
correlaliuns. LaLei, X-ia^^ ubsei vaLiuns uf giuups uf galax- 
ies have strengthened the case for a non-gravitational en- 
tropy injection in these systems, in particular by the ob- 
servation of the so-called "entropy floor" (Ponman et al. 



et al. (1997, 199S, 199^, hereafter CMT) have introduced 



199£). These authors showed that the central entropy level 
in groups (outside the cooling-flow radius) is higher than 
the level the sole cosmological shocks can provide, while 
being well in accord with adiabatic numerical simulations 
[e.g., including no dissipative physics) in clusters. This en- 
tropy injection could be due to supernovae explosions (the 
so-called feedback from star formation) or other sources as 



a simple model, where the bending is provided by the dif- 
ferential strength of a shock occuring at the virial radius 
in a preheated gas. Their description relies on the physi- 
cal modelling of the shock interface, using Hugoniot rela- 
tions, and reproduces naturally the central observed en- 
tropy floor. While this model highlights the important role 
of shocks and entropy floor in the formation of clusters and 
groups, no analytic description of the X-ray observed re- 
lations has yet been found to encompass both groups and 
clusters. The aim of this paper is to present such a model 
and compare it to observations. 
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However, despite the success of these semi-analytic 
schemes when confronted to observations, difficult prob- 
lems remain, as for example the fact that a reheating by 
supernovae (hereafter SN) explosions requires an incredi- 
bly high efficiency of the transfer of energy from SN rem- 
nants to the intergalactic medium (Valageas & Silk, 1999; 
Bower et al., 200"o| ). In particular, the cooling of the SN 
remnants will obviously decrease this transfer efhciency, 
and must be ignored in these models. While a combination 
of SN and quasars (hereafter QSOs) reheating is probable 
and would alleviate this problem, Bryan (200C) has pro- 
posed that differential galaxy formation between groups 
and clusters (which would lower the central entropy in 
groups and allow higher entropy gas to flow into the cen- 
ter) can explain the curvature of the Lx — T relation and 
the entropy floor. However, it is well known that with- 
out preheating, most of the available gas in the universe 
would have formed stars today, which is not observed (this 
the so-called "overcooling problem" , see Blanchard et al 



1992; Balogh et al., 2001 ). This much- needed feedback in 
galaxy formation would certainly have an impact on the 
formation of clusters and groups. Valageas & Silk (1999) 
have indeed shown that galaxy and cluster feedback differ- 
ing requirements are likely to provide a tight constraint on 
the amount of preheating in the universe. But, even in the 
preheated scenarios, there is actually no consensus on the 
entropy injection epoch. While most of the studies have 
focussed on "external preheating" models (where the en- 
tropy injection occurs before the formation of groups and 
clusters e.g., before the gas is compressed by shocks) be- 
cause low density gas requires much less energy than high 
density one to be raised to a given entropy level (and also 
because the star formation rate seems to peak rather early 
in the universe evolution) , [Loewenstein (2000 ) , from a s c- 



ries of approximat e static hot gas models and Brighcnti 
& Mat hews (2001 ), from ID numerical simulation with 
cooling, mass dropout and star formation feedback, have 
argued that most of the heating occured during or after 
the assembly of the group or cluster gas. The efficiency 
problem of SN explosion is still present, but can be allevi- 
ated by an initial mass function flatter than the Salpeter 
one ( Brighcnti fc Mathews, 2001 ). These are the so-called 
"internal preheating" models. 

In this paper, we focus explicitly on an external 
preheating model. Internal preheating models require 
spatially-dependent and time-dependant star formation 
rates to provide the amount of injected entropy, while ex- 
ternal preheating models only require the level of entropy 
before the gas falls in a cluster or a group. This simplifica- 
tion allows us to obtain a completely analytic model that 
describes both groups and clusters. The physical approach 
we follow consists of deriving scaling relations, e.g. Lx—T, 
y — T, y — Lx by linking these observed quantities to the 
specific entropy profile in the system. For this purpose we 
first derive an analytic expression for the entropic profile 
normalisation at the virial radius of clusters and groups 
by discussing shocks, as in CMT. Using this relation we 
then derive analytically the relevant scaling relations and 



compare them to data. Following Tozzi & Norman (2001), 
we highlight the key role of entropy and argue that it in- 
deed constitutes the best observable in clusters, allowing 
to derive analytical expressions for standard X-ray and 
Sunyaev-Zeldovich (hereafter SZ) correlations. This leads 
us to reinterpret these relations in terms of global entropy 
content in a well physically motivated manner. 

The plan of the paper is as follows : in section |^, we 
present the shock model and derive an expression for the 
normalisation of the entropy profile at the virial radius. 
This expression, containing a free normalisation, is fitted 
to observations. In section ||, we relate the X-ray luminos- 
ity to the entropy profile and derive an analytic expression 
for the shape of the Lx—T relation, which is subsequently 
compared to observations using only the fitted parameter 
of the last section. The same method is followed in sec- 
tions ^, ^ and 1^ to obtain analytic expressions for the 
MicM - T, fgas - T, yo - T and j/q - L^ correlations 
(where Mjcm is the gas mass, fgas is the gas fraction and 
yo is the central Compton parameter). In each of these 
parts, we provide a detailed comparison with previous 
work. Section ^ discusses a reinterpretation of these corre- 
lations in terms of the total entropy content in groups and 
clusters, tests it using a semi-analytic model and discusses 
the main hypotheses we make. Moreover, a comparison 
of our results for the entropy floor with other theoretical 
models is made. Section ^ summarizes the present work 
and concludes. We also compute the numerical values of 
the normalisations of the correlations found and summa- 
rize all these results in Appendix allowing to use them 
easily in another context. Appendix ^ computes the in- 
falling velocity as a function of the mean mass of the sys- 
tem and compares to hydrodynamic simulations. Finally, 
Appendix introducing a simple model for the entropy 
profile of groups, computes the shape factors appearing in 
the different normalisations. 

Unless otherwise stated, we use Hq = 
100 /i2/3 km s^^ Mpc~^ with ^12/3 = 2/3, which gives 
approximatly Hq — 67km s^^ Mpc^^. We rescale all the 
data used in this paper to this value. 



2. Entropy and shocks in clusters and groups 

The study of the equilibrium state of X-ray clusters and 
groups of galaxies requires the physical description of 
the complex interplay between baryonic and non-baryonic 
dark matter. To tackle this problem analytically requires 
several simplifying hypotheses. The following assumptions 
are reasonable and will be thus assumed henceforth (see 
^arazin, 1988 , for a review and a discussion of the main 



assumptions). 

First, the hot and diluted plasma of ions and electrons 
constituting the intracluster medium (hereafter ICM) is 
treated as a single species perfect gas satisfying 



P = 



pkT 
firUp 



(1) 
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where P, p, T denote respectively the pressure, density 
and temperature of the baryonic gas, while k denotes the 
Boltzman constant. Second, neglecting the ICM mass with 
regards to the DM mass, and assuming stationarity (no 
gravitional potential variation on time scales smaller than 
the hydrodynamic time scale, e.g. no recent mergers), the 
ICM is assumed to be in hydrostatic equilibrium (hereafter 
HSE) in the dark matter gravitational potential (j>. We 
can thus write the continuity and the Euler equation as 
follows : 



obvious in clusters from X-ray observations for example) 
and, locally only, the system is in thermodynamical equi- 
librium. Thus, the state functions have the same depen- 
dency on state variables as in global equilibrium, but vary 
now from point to point as do the latter. Note that this 
hypothesis is likely to break down at smaller scales than 
the ones of interest to us, due to turbulence or magnetic ef- 
fects. This hypothesis is crucial if we want to use the usual 
analytic expression of the specific entropy of the ICM, s, 
namely 



VP 







(2) 

pV0 (3) 

where v is the gas velocity and the gravitational po- 
tential. Equations (||), (||) and (|l|), together with the en- 
ergy conservation equation and Poisson equation (the use 
of which can be avoided if an approximate analytic ex- 
pression for (j) is chosen, e.g., neglecting the gravitational 
contribution of the baryons) form a closed set of equa- 
tions, whose resolution provides the radial profiles of all 
the quantities of interest in the gas. 

To solve this set, boundary conditions are required. In 
the present paper, we are interested in the variations of 
these boundary conditions with the total mass of the sys- 
tem (or equivalently with its mean temperature). Indeed, 
the previous equations will describe the equilibrium of 
baryonic gas in the underlying potential, whatever the to- 
tal ma ss of the system. We argue and show in the following 



S/ Cy — In 



kT 



(4) 



where 7 is the polytropic index of the gas. The thermo- 
dynamical state of the gas is fully known as soon as the 
entropy profile of the gas is known. Since we want to high- 
light the key role of entropy as "an observable" we choose 
to describe the physical state of clusters in terms of en- 
tropy, or more appropriately and equivalently in terms of 
the adiahat, defined as 



K = (/imp) 



7-1 , 



kT 



{kT)n 



-2/3 



(5) 



where n = p/{iimp). In the following, T will denote the 
temperature in keV (i.e., we will replace everywhere kT 
by T). Note that from now on the polytropic 7 index will 
be fixed to its standard value of 5/3. The definition of the 
adiabat differs from the widely used definition of Balogh 



that c hanges in the boundary conditions, when explicitly 
written as a function of the mean temperature, allow us 
to derive the change in global correlations, such as the 
Lx — T relation, on a scale going from clusters to groups 
of galaxies. For this purpose, we use the fact that the flow 
of gas on a cluste r becomes supersonic in the vicinity of the 
virial radius r„ ( Tcyssier et al., 1997 ) . We will then show 
that a standard modeling of the resulting shock (following 
CMT) leads to the appropriate boundary conditions. We 
will consider hereafter that the boundary conditions take 
place after the last major merger of the system, as appro- 
priate in a hierarchical picture of structure formation. 

The last hypothesis we will use is that within the virial 
radius, the ICM is isothermal. This hypothesis will not be 
used in the derivation of the boundary conditions (sec. 2.1) 
but in the derivation of the global correlations in groups 
and clusters (sees ^ and ^) . Both observations [Markeviteh 
et al. ( 1998| ) and simulations (Teyssier, 2001) show that 
this approximation is good into a factor two and is thus 
sufficient for our purpose. 

Finally, let us recall an hypothesis necessary underly- 
ing every work published using the entropy of the ICM, 
but rarely stated : the local thermodynamic equilibrium 
hypothesis (hereafter LTE). By definition, in a state of 
global thermodynamic equilibrium, temperature and pres- 
sure (and thus density) are constant throughout the sys- 
tem, and the state functions of the system (e.g. entropy) 
have a determined dependency on the state variables cho- 
sen (and so a fixed value). In the LTE however, temper- 
ature and pressure can vary from point to point (as is 



et al. (1999) (by a constant factor {p.mp)'^ ^) but allows 



a direct comparison to observations. 



2.1. A general shock model 

We first aim at deriving adiabat boundary conditions 
based on the Rankine-Hugoniot relations, used in the same 
form as CMT. 



2.1.1. The Rankine-Hugoniot relations 

Let the infalling gas velocity reach the sound speed at 
r — ry and consider subsequently that a shock form at this 
radius. Let the temperature, the density and the velocity 
of the infalling gas before the shock be respectively Ti, 
ni and Vi, and the post-shock temperature, density and 
velocity be T2, 71.2 and W2- If the shock efficiency is very 
high, i.e. the post-shock velocity V2 is null in the rest- 
frame of the cluster center of mass, the Rankine-Hugoniot 
relations yield (CMT) : 



fcT2 = ^F(e), 
with 



F(6) = 



7_ 3 

10 ^ 20 (1 + vrny 



and 



15 kTi 
4 /i nip vl 



(6) 



(7) 



(8) 
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At the same time, the ratio between post-shock (rt2) and 
pre-shock (ni) density is given by : 



4 1- 



To 



Ti 

T2 



1/2 



(9) 



Note that since the ICM within is assumed to be 
isothermal, its temperature is equal to the temperature 
T2 after the shock. 

Using the formulas above, we will now derive a gen- 
eral expression for the adiabat jump at the shock, which 
will depend on the system mean temperature. This virial 
adiabat normalisation will provide us with the change in 
the boundary conditions necessary to describe groups and 
clusters in a unique analytic model. 

2.2. The adiabat after the shock 

To find a general expression for the post-shock adiabat, 
we first express the ratio T1/T2 as a function of e : 



Ti 
T2 



5F(e)- 

Introducing this expression in equation 

4 e X 



(10) 



one has 



712 — 2r7,i 



+ 



1 - 



5 Fie)' 



5F{e)J 5F{e) 



1/2 



(11) 



Using then equations (^), (||) and (|Tl|), the post-shock 
adiabat K2 is easily expressed as a function of the pre- 
shock adiabat Ki and e in the following manner : 



K2 



kT2 
2/3 



(12) 



where 



Hie) = 



F{e) 



(( 



1 - 



5F{e) 



+ 



5F(e) 



5 Fie) 



1/2 



-2/3 



(13) 



Even if this last expression looks complicated, it can be 
straightforwardly expanded in a Laurent serie in the vicin- 
ity of e = 0. We obtain : 



Hie) =. II 2-/3 (1 + 1, r 



Oie) 



(14) 



Checking the validy of this expansion in the range 
< e < 1 we find an agreement better than 4% (fig- 
ure |l|) which is enough for our purpose (the error of the 
approximation reaches a constant 15% level at large e). 
Thus we will neglect higher order corrections. 



2 3 
f 2 




Fig. 1. Top : plot of the function iJ(e) (equation ^ solid 
line), together with the two firts order of its expansion 
(equation |lj, dashed line) . Both functions are hardly dis- 
tinct. Bottom : Percentage difference between H and its 
expansion. The agreement is better than 4 % in the rele- 
vant e range (and always better than 15 % for all values 
of e.) 

As a consequence, for the all range from groups to rich 
clusters, one can write to a very good approximation : 



17 i/c, / 10 1 
16 V 17 e 



(15) 



Note that this expression that fits so well the general adi- 
abat expression is not really a surprise. Its physical signif- 
icance is indeed straightforward. To show this, let us take 
two limiting behaviour, namely in the cluster and in the 
group mass range. 

Since rich clusters of galaxies accrete mainly small 
clumps and diffuse gas, the cold inflow or strong shock 
limit is appropriate. Consequently it is usually argued that 
the entropy of rich clusters is dominated by shock created 
entropy. Indeed, for this approximation we have T2 ^ Ti 
and e ^ 1, and so (from eq. |l5|) 



5 2-1/3 
K2^- i^i. 



(16) 



On the other hand, for groups of galaxies, the weak 
shock or adiabatic infall approximation is more appropri- 
ate. Indeed since the infall speed tends to be lower, and 
since groups accrete preheated gas, the weak shock limit 
(whose limit is the true adiabatic infall) is appropriate 
(CMT)[[ Indeed, taking the limit e ^ 1, eq. (^ yields 



27 
16 



2-1/3 



(17) 



This means that the entropy content of groups is domi- 
nated by the "adiabatically accreted" gas, hence the ex- 

^ Note that this hypothesis of adiabatic infall has been shown 
to allow to reproduce the Lx — T relation in the group range 
by [Balogh et al. (1999|) 
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istence of a so-called entropy floor (the last expression of 
K2 being independent of e). 

The analytic expression in eq. ( |l5| ) thus provides a 
well physically motivated expression for the competition 
between shocks and entropy floor that rules the adiabat 
virial normalisation from groups to clusters. 

2.3. Linking groups and clusters 

2.3.1. Defining a general analytic K — T relation 

We now want to make use of the previously derived ex- 
pression ( |l5|) to predict a general relation for the K — T 
relation for groups and clusters. 

First, making use of eq. (0) and eq. (fsh we rewrite 



if2 



-2-1/3 /f, 
16 ^ 



1 



8 /i TOp Vi 

51 kTi 



(18) 



Second, we need to express the infall velocity Vi as a 
function of the mean temperature of the system. This is 
done in Appendix where we show that using the virial 
theorem, the M—T relation, together with the assumption 
of adiabatic cold inflow, we can write {jpi being a caracter- 
istic virial radius gravitational potential, T* a caracteristic 



temperature and 77 a dimensionless constant, see eq. B.5 
B.ll| and |B.13D : 



(19) 



Note that from now on we will assume notations of 
Appendix 

This leads naturally to: 



^ ^ 10 



16 nmp(j>i 

1 77 — 

51 ' kTi 



that we rewrite 



1 



-) 



where we have defined^ 
17, 



Ko 
and 



16 



2-1/3 K, 



51 kTi 



(20) 



(21) 



(22) 



(23) 



16 ?] iinip 01 

The scaling relation ( |lj[ ) is a key result of our paper. 
It states how the entropy jumps due to the shock changes 
with the mean temperature of the system. Written this 
way, the last formula has two free parameters : a normal- 
isation parameter Kq ~ Ki^ i.e. the adiabat before the 
shock, and the temperature Tq that marks the transition 
from the adiabatic inflow to the strong shock regime. 

Deriving the value of Ki requires a whole semi-analytic 
scheme, which would require a model of the entropy evo- 
lution of the intergalactic medium (hereafter IGM) as the 

1/3 



universe evolves (see e.g. Valageas & Silk, 1999). This is 
clearly beyond the scope of the present paper, and we 
will thus constrain Ki by comparison with observations 
(which only means that the normalisation of our model is 
left free). However, Kq and Tq are not independent param- 
eters and can be related by simple physical considerations 
as will be shown in the next section. 

Once this is done, the scaling relation ( pl| ) is left with 
only one normalisation parameter Kq- 

2.3.2. Relating Kq and Tq 

First, to explicit the relation between Kq and To, we need 
to specify the value of both r] and Ti in equation (p3|). 
Using the results of Miniati et al. (200C) for the infall 
velocity in hydrodynamic simulations (see eq. B.15 and 



Appendix for notations), a value for rj can be derived by 
equating eq. (B.15) and eq. ( pj| ) : 



which, when introduced in equation (|23|), gives : 



Tn = 



51 kTi 



■T, 



(24) 



(25) 



8 umpV^Q 

where Ts and Vsq are known (appendix ^). 

The relation between Kq and Tq will appear when ex- 
pliciting Ti. Indeed we can write it as a function of the 
pre-shock adiabat Ki : 



Ti 



2/3 



(26) 



The overdensity that eventually collapses and reaches 
a density pi at has decoupled from the Hubble ex- 
pansion at a turn-around redshift zta whe re its density 
was p ta and its radius was rta- Following Balogh et al 



|(1999|) , we assume first that the IGM was preheated (non 
gravitationally) before zta and its adiabat raised to the 
value Ki . Second, we also assume that between rta and the 
shock radius r^ itself, the gas has an isentropic behaviour. 
To fix the value of the overdensity (and thus the density) 
^ta — [pta — p)/p at tbe turn-around, we will consider the 
simplest collapse model, i.e. the spherical collapse model 
in a given cosmology. 

By definition of the baryonic fraction fs and the over- 
density Sta we can write : 



Pta = (1 + Sta) Pc fh 



(27) 



where pc is the critical density. Inserting this last equation 
in eq. (^), one obtains : 



To 



51ifl {{l + 5ta)Pc.fBf''' 



8 



(/ITOp 



(28) 



Remark that 17/16 x 2" 



: 0.843, so that Ko ~ 0.843 K 



which relates Kq and To via eq. (|22|). 

Consequently, to fix the numerical value of To we need 
to specify a value of Ki. We choose to assume the value 
of the observed entropy floor visible in groups of galaxies. 
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A value of Ki = 120kcVcni^ is inferred from the data 



(Lloyd-Davies et al., 200C, see the next section). 

Then we assume a flat = 1 universe with 

Hq — 50kms"^Mpc~^ and the baryonic fraction as 
fixed by Big-Bang Nucleosynthesis standard, i.e. fs — 



0.015 h-^Qf^/Vlm ( Olive et al., 2000| ). For a spherical collapse 



model 5ta — 97r'^/16. Note however that in the widely used 
hierarchical structure formation paradigm, this monolithic 
formation model is too simplistic. Indeed, since the ma- 
jority of the infall of group of galaxies occurs through fil- 
aments, an overdensity of Sta ^ 10 should be more appro- 
priate. Considering then the values of Vgo and spec- 
ified in the Appendix, we find that Tq = 1.86 keV for 
= 97rVl6 and Tq = 2.62 keV for 6ta = 10. This is 
quite insensitive to the adopted cosmology, since for an 
= 0.3 universe, one finds Tq = 1.71 keV and 2.33 keV 
respectively. It is thus obvious that we can take for Tq an 
intermediate value between these two extremes. 

In the following, we will thus consider that Tq is fixed 
at 2 keV, and that the only free parameter of the model 
will be the normalisation factor Kq. Moreover, choosing 
a value of Tq between 1.5 and 2.5 keV does not alter the 
qualityi of the fita we present in the next scotiona. 



Looking aL eq. Tq has a shuple physical luLui'pi'c- 
tation : it is the system mean temperature where the ef- 
fect of an entropy floor begins to change the adiabat jump 
during the virial shock. For systems with T < Tq, the pre- 
shock entropy is high enough to reduce the shock and a 
quasi-adiabatic inflow takes places (but the shocks are not 
completely suppressed, see section 7_A). For T > Tq, the 
shock adiabat jump brings the post-shock adiabat to val- 
ues much higher that the entropy floor, the latter having 
less and less incidence as T increases. Thus we find that 
this transition occurs around 2 keV. 



2.3.3. Comparing with observations 

X-ray observations of groups and clusters of galaxies have 
brought a great amount of information on the gas equi- 
librium in these systems. Unfortunately, due mostly to 
the high background of the instruments and the gas den- 
sity decline with radius, there are not yet entropy obser- 
vations available up to the virial radius, that would al- 
low to test directly our model (and more generally the 
scheme invented by CMT). The new generation of satel- 
lites (Chandra and XMM-Newton) is expected to bring 
new observations near the virial radius, allowing to test 
different formation scenarios of clusters and groups (see 
e.g., [Tozzi et al., 2000|). 



Nevertheless, some data on the central adiabat value 



of the gas in groups and clusters are available. Ponman 



et al. (1999) were the first to show that groups have a 



central adiabat higher than the one they would have if it 



had been imposed by shocks only. Recently, Lloyd-Davies 



et al. (20001) have refined these observations, taking into 
account the temperature gradients, while increasing the 
statistics. This is the so-called observed "entropy floor". 




T (keV) 



Fig. 2. Post-shock adiabat K2 as a function of the mean 
temperature of the system. The solid line is obtained 
with equation (pT|), with Ki = 120keVcm^ {i.e., Kg = 
102keVcm^). Note that only the normalisation is a free 
parameter, as the reference temperature Tq = 2 keV is ob- 
tained from the value of Ki . The dashed line is a X cx T fit 



to the systems above 4keV found by Helsdon & Ponman 
(2OOOD. The similarity of the shock model adiabat nor- 



malisation (computed at the virial radius) with the data 
taken at the center {R = 0.1 r^,) allows us to postulate an 
approximate self-similarity in the specific entropy profile 
in the Sec. 0. 



While cool systems have a constant central adiabat, hot- 
ter clusters follow naturally the results of numerical sim- 
ulations without a preheating phase (which predict that 
K (X T). Note that our model predicts exactly the same 
behaviour, exhibited in the post-shock adiabat analytic 
expression of eq. (^l|). However, we compute the adiabat 
near the virial radius, while the data come from the cen- 
tral parts of the systems. Nevertheless, observations and 
theory can be compared if we assume that the adiabat 
difference between two different temperature systems is 
conserved when going from the virial radius to the central 
part. This implies that the physical mechanism responsi- 
ble for the post-shock adiabat value (competition between 
shocks and an entropy fioor) is the same as the one re- 
sponsible for the central adiabat value. This is indeed a 
reasonnable assumption, since 3D hydrodynamical simu- 
lations of the formation of clusters show that, after the last 
major merger, a quasi-spherical shock forms and expands, 
that leaves behind an increasing adiabat profile (see e.g., 
Evrard, 1990; Frenk et al., 1999). The same physics is 



thus at work in the center and in the outer parts of clus- 
ters, allowing the physical description of the two regions 
to coincide. Indeed, the only difference between these two 
regions is that the universe will have evolved when the 
shock reaches the present-day virial radius, as compared 
to the time when the central entropy is set. We thus don't 
expect the observational value of Ki (found by normal- 
isation to the observations) to be representative of the 
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present adiabat of the universe, but it should instead be 
representative of the IGM at the time when the smaller 
mass systems formed. If virial radius observations of the 
adiabat were available for the same systems, we expect 
then that the same analytic expression as in eq. (|^) would 
be a good description of the data, with a different value of 
Ki, representative of the present-day adiabat of the IGM. 

F igure |^ shows the compari son of eq. ( pl| ) with data 
from Lloyd-Davies et al. (200(j| ). Only the normalisation 
was fitted here, with a value of Ki = 120keVcm^ (giv- 
ing Kq = 102keVcm^), while wc have taken Tq — 2keV, 
derived with the same value of Ki in the last section. 
Replacing eq. (28) in eq. ( |2l| ) instead, and fitting for Ki 
gives an undistinguishible analytic curve. The dashed line 
is the result of a if oc T fit to the systems above 4keV 
found by [Helsdon fc Ponman (2000 ). This temperature 
dependence of the adiabat is ta ken from numerical simu- 
lations (see Ponman et al., 19991 ) . The dot-dashed line is a 



constant ga s entropy fitted to the four lowest temperature 
systems by Lloyd-Davies et al. (2000 ) and has a value of 
139keVcm^. The agreement between the analytic formula 
and the data is remarkable. Indeed, using the fitted value 
of Ki, the computed Tq value agrees very well with the 
temperature of the intersection between these two limit- 
ing behaviours (see last section). This shows that eq. 
is not a mere analytic formula (which could have been 
infered from the observations), but succeeds in capturing 
both ingredients which fix the central entropy value : the 
entropy floor (dot-dashed line in the groups mass range) 
and the shocks (dashed line, in the cluster mass range). 

By adjusting the normalisation of our analytic model 
using the data we derive a value for Ki . Then making use 
of this relation we get directly for any system at a given 
temperature T (either a group or a cluster) the normali- 
sation of the adiabat profile. It is now important to check 
wether we are able to use this relation to derive realistic 
Lx — T and y ~T relations for both groups and clusters. 



3. From the adiabat to the Lx — T relation 

The total X-ray luminosity of a local cluster (group) of 
galaxies is expressed as : 



Lx= I nlA,{AE,Te)dY, 
Jv 



(29) 



where Ue is the electronic density and the integration 
is carried out over the whole volume V of the system. 
Ae(A_E, Te) denotes the local X-ray emissivity of the clus- 
ter gas within a given energy band AE. For example, in 
the case of pure thermal Brehmsstrahlung emission ob- 
served with contemporary wide-band satellite, one has 



A, 



^1/2 



3.1. An analytic prediction for the Lx — T relation 

We can now express equation (|29| ) as a function of K and 
thus link the X-ray luminosity with the temperature and 
the adiabat : 



X3 



Ae(A£;, T) dV 



(30) 



Since we assume the cluster to be isothermal and since 
in general the emissivity can be written as a power-law 
function of the temperature with index a, one has : 

dV 

K3' 



Lx 



An 



(31) 



Aq being the normalisation of the emissivity. 



Equation ( pl| ) gives us a direct link between Lx , T and 
K. Was the expression of iiT as a function of T known, it 
would enable us to predict the shape of the Lx — T rela- 
tion. However, the result of the previous section (^l|) gives 
us the change in the adiabat normalisation as a function 
of temperature, but since the integral in eq. (^) is com- 
puted over all the system volume, in principle we need the 
adiabat profile to compute it. However, if we make the fur- 
ther assumption that the adiabat profiles are self-similar 
in temperature, i.e. the temperature enters the analytic 
expression of K[r) by its normalisation only, then we can 
derive an expected Lx — T relation. This assumption is 
justified by the fact that our theoretical expression for the 
adiabat normalisation (equation (pl|), which is computed 
at the virial radius, provides a very good fit to the cen- 
tral entropy data in groups and clusters. This simple fact 
ensures that the difference in normalisation between two 
clusters of different temperature (which is due to shocks 
and preheating) is leaved approximately unchanged from 
the virial radius to the center. This important side aspect 
of our work will be discussed in Sec. 0. We thus write : 



K{v,T2)^K2{T2)f{v) 



(32) 



where K2{T2) is given by equation ( ^l| ) and /(r) is a func- 
tion of the radius, that is independent of the mean system 
temperature. We have thus : 



^=Aoi^2(T)- 



dV 

7^' 



(33) 



Even if / does not depend on T, the integral on the 
right-hand side of the last equation depends on the total 
volume of the cluster, and thus on its virial radius, which 
depends also on the mean temperature. Let us assume that 
the cluster is spherically symmetric (which is reasonable 
if the system has sufficient time to relax). We can then 
write : 



dV 



Att / r^f-^{r)dr 
Jo 

Airrl / x^f^^{x)dx 
Jo 

„3 



4:Trr^ Qx 



(34) 

(35) 
(36) 
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that we can compute a value for Qx by considering 
aa reasonable model for the adiabat profile, relying on 
isothermality and a (3— model for the gas density profile 



(Cavaliere & Fusco-Femiano, 1976). The computed value 
is Qx = 0.115. 

To validate this relation, we plot on figure ||the Lx — 
T relation of equation (|3^) together with data taken for 



groups (f rom Helsdon fc Ponman, 2000) and clust ers of 
galaxies ( |Markevitch, 1998| ; [Arnaud fc Evrard, 19991) . The 



T (keV) 



Fig. 3. Lx — T re lation versus obs ervations from Helsdon 
fc PonSan (ZUUu|j tor gr oups and LVlarK:evitcli (liJiJS ) and 



quality of the fit is obviously excelle nt. No te that we adopt 
for Tq the value computed in Sec. 2.3.2 and used Qx = 
0.115. There is thus no free parameter in the solid curve 
computation, since has been already fixed. 

Note that the hotter groups seem to be overluminous 
when compared to the cluster data as well as to the ana- 
lytic relation. This is due to the fact that, in the cluster ob- 



servations, the central cooling-flows were cut ( Vlarkevitch 



Arnaud k. iivrard (lyyb ) tor clusters, ihc solid fine is the 
analytijc tormula m cq. ([?y|). Jq was nxcd at the value 

UtE 



1998) or the sample was chosen explicitly to contain 



computed m Sec. 'i.A.,\ i.e., ip = ^ ke V , and Ljx = U.iiQ, 



as computed in Appendix ^ 

where we set x ~ r/r^ and define the shape factor Qx- 
Using then the — M and M ~ T relation of equa- 
tion dB.lCl ) and ( |B.ll| ) before using our K — T relation 
(^), we obtain when replacing in equation (|3^) : 



only clusters with known weak cooling-flows ( Arnaud fc 
Evrard, 199£). On the contr ary, the groups luminosity 







1 + 


(^)l 







-3 



where the normalisation Lq is given by : 
3 AL 



Lo = Ao 



3/2 



Ko'Qx, 



(37) 



(38) 



Let us first comment on the shape of the Lx — T relation. 
While for the group range, where a ~ and T Tq, we 
have : 



Lx oc r9/2 cx ■^ 



(39) 



for very hot clusters (for which T ^ Tq and a = 1/2), we 
recover the so-called self-similar behaviour : 



Lx oc r9/2+l/2-3 ^ y2 



(40) 



The 10 keV cluster range will be intermediate between 
these two behaviours. We now turn to a comparison with 
observations. 

3.2. Comparing with Lx — T observations 

The normalisation of the local Lx ~ T relation (eq. ^) de- 
pends on the value of the IGM adiabat Kq and the s hape 



in Helsdon fc Ponman (200C ) are not corrected for the 
cooling-flow (which explains that the hotter groups, which 
are expected to have the greatest fraction of their luminos- 
ity coming from their cooling-flow because of their higher 
central density, show the greatest departure). Our model 
does not take into account a possible cooling-flow compo- 
nent, which explains that it fits very well the cluster data, 
while lying near the lower envelope of the hotter groups. It 
fits nevertheless very well the lower temperature groupg|^ 

Both the cluster slope and the group steeper slope 
come out naturally of our relation. This agreement is re- 
markable and is the first to our knowledge based on an 
analytic discussion to encompass both groups and clusters. 

This good agreements is a strong sign in favor of the 
consistency of our approach. It is thus natural to try to 
apply an analogous method to the correlation between the 
mass of gas and the mean temperature, as well as to the 
SZ effect, and compare cluster observations with our ana- 
lytic predictions. Unfortunately, gas mass observations are 
sparse and heterogenous for groups, while no SZ effect in 
a group has yet been observed. We can nevertheless make 
predictions for the Mjcm — T and y — T relation in this 
range, which could be validated by undoubtly forthcoming 
SZ measurements. 



4. From the adiabat to the ICM gas mass and gas 
fraction 



4.1. Predicting the Mjcm 
to observations 



T relation and comparing 



factor Qx- The former has been obtained in Sec. |2.3.3 
{Kq = 102keVcm^) by comparing our analytic formula 
for the adiabat jump to observations. The latter is more 
difficult to obtain since it will depend on the exact adia- 
bat profile of clusters and groups. We show in Appendix |c| 



Another tight correlation found in reducing ICM data of 
clusters is the link between the ICM mass and the mean 
temperature. For a perfectly self-similar model, one ex- 
pects theoretically Mjcm oc T'^/'^. In fact, several groups 



^ Fo r more discussion on this point, see Brighenti < 
1(200 1[ ), Sec. 2.1. 
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have reported recently a steepening of this relation, which 
could be due to entropy injection. Vikhlinin et al. (1999) 
have studied the outer regions of a ROSAT sample of clus- 
ters, and found that the above correlation could be writ- 
ten MicM oc yi 'i'iio.is^ steeper but nevertheless close 
to the self-similar value. On the other hand, Mohr et al 
(1999|)]have fitted double /3-models to another large sam- 
ple of clusters and found that Mjcm oc t^-^s,±o.18 ^ much 
steeper than the self-similar model. It is not clear if this 
discrepancy (at only 1 a level) is real and if so, what are 
the reasons behind it. 

It is easy to compute this correlation in our framework 
and worth comparing our predictions with observations. 
The ICM mass of gas can be written (assuming spherical 
symmetry) : 




T (keV) 



M 



ICM 



Ajl 



(41) Fig. 



et al 



4. Mjc M — T correlation. The data are from Mohi 
. ( 1999| ) . The dot-dashed line is the best-fit linear cor- 



where Pgas{r) = fj,mpn{r) is the mass density of plasma. 
Rescaling with the virial radius and replacing the gas num- 
ber density by a combination of temperature and adiabat, 
one finds : 



MjcM — Mj^ T 



1 + ^ 
To 



1 -3/2 



with 



M 



ICM 



/i Trip 



3Af, 



3/2 



-3/2, 



!M, 



where we defined naturally 



;m 



^ dx 

7^ 



(42) 



(43) 



(44) 



The value of this shape factor, which determines the over- 
all normalisation, is also evaluated in Appendix |^. The 
predicted value using eq. (^ ) is Qm = 0.1915, while 
the best-fit normalisation is ^ 40 % lower. We argue in 
Appendix ^ that the observed surface brightness in most 
of the clusters does not reach rsoo (the virial radius) and 
thus needs an risky extrapolation of the data. Moreover, 
the observational derivation of ta is also dependent on 
several assumptions (in particular on the assumption of a 
given mass profile) and it is possible that the gas masses 
were in fact computed at A 7^ 500. Indeed, lowering the 
the upper boundary of the integral in eq. (^) to 
{i.e., assuming rsoo is underestimated by only 20%) gives 
Qm — 0.1136, which is in perfect agreement with the data. 
We use this value in figure ^ and in the following section. 

The slope of the correlation also steepens from very 
hot clusters (Mjcm oc T^/^, analogous to the self-similar 
correlation) to groups {Mjcm oc T^). The temperature 
range probed by Mohr et al. (1999| ) is intermediate be- 
tween these two behaviours. 

The figure | shows the data from |Mohr et al. (1999|) , 
as well as their best-fit linear correlation Mjcm oc T^-'^^ 
(dot-dashed line). The sohd line is our predicted relation 
using the computed value Qm — 0.1136. The predicted 



relation Mjcm oc T^-^^ found by these authors, while the 
solid line is the predicted correlation (eq. p^) using the 
computed value Qm = 0.1136 (see Appendix |C|) . The cur- 
vature of our prediction is clearly visible when compared 
to the best-fit power-law. 

slope matches perfectly the data. The steepening of the 
relation compared to the self-similar one is thus, in our 
particular model, simply a consequence of the differential 
shock strength and the entropy floor. 



4.2. Comparing the subsequent fgas 
observations 



T relation with 



The gas fraction is very important in that it is supposed 
to give a direct lower limit on the universe baryon frac- 



tion, if clusters are a fair sample of the universe (White 
[st al., 1993| ). Such measurements have proven to be tight 
constraints to the total mass density of the universe. 

Using eq. (^2|) and the observed M — T relation of 
eq. ( B.ll ), i.e. the observed scaling between total mass 
and temperature, it is straightforward to compute the gas 
fraction variation with temperature : 



f - fO 7-3/2 

J gas — J qas 



8 with 

fgas ^ ^^mp 



1 



T 



-3/2 



-3/2 



Qm 



(45) 



(46) 



The value of Qm is the same as in the last section {Qm 
0.1136) and no renormalisation is performed. 

Eq. (^H) is shown in figure |5|, together with observa- 



tions. The data points are from p\/Iohr et al. (1999 ). These 
data were fitted by the authors with a power-law, namely 
fgas = (0.207±0.011)(T/6keV)°-34±0-22. The shaded area, 
limited in abscissa by the minimum and maximum tem- 
perature of the sample, represents the area covered by 
similar power-laws increased and decreased by 1 (both 
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T (keV) 



Fig. 5. fgas — T relation. Data are taken from Mohr et al 
(1999|}]as well as the shaded area which corresponds to 
the domain of validity of their fitted power law, ±1ct. The 
solid line corresponds to our prediction. Both the trend 
and the normalisation are in good agreements. 

in normalisation and in slope) . It is obvious that the data 
are dispersed, but the trend found by Mohr ct al. (1999 ) 
is very well recovered by eq. (^5|), as well as the normali- 
sation. 

Looking at eq. (45), one sees that for Tq <^ T, i.e. for 
rich clusters, the gas fraction becomes independent of the 
temperature. This regime is only reached asymptotically 
in the figure. On the other end the slope steepens from 
clusters to groups, since we have : 



fgas oc fg^g Tq ' 

cx M 



constant for rich clusters. 



for groups, 



(47) 



where M is the total mass of the group. It is interesting 
to note that the fact that the gas fraction is proportional 
to the total mass in the groups regime was first claimed 
by Balogh et al. (1999 ). We thus recover analytically the 
adiabatic behaviour modelled by these authors. However, 
their model breaks down at high temperature, when the 
shocks become important, while ours runs smoothly from 
an adiabatic infall to a shock-dominated regime. 

5. From the adiabat to the SZ temperature 
decrement 

The observations against which our model was validated 
above were only in the X-ray wavelengths. We are now 
turning to the radio and sub- millimeter bandwidth, which 
also probes the hot gas in clusters (via the Sunyaev- 
Zeldovich effect) and offers independent observations, with 
different possible systematics and errors. This allows us to 
independently validate the model and make some predic- 
tions about the group behaviour in this band. 

Inverse- Compton scattering of incoming CMB photons 
on ICM thermal electrons cause a well observed change 
in the spectral dependence of the CMB ( Rephaeli, 1995 ; 



Birkinshaw, 1999). The amplitude of the CMB tempera- 
ture brightness variation can be written : 



AT 

TcAIB 



= h{x) 



kT 

meC^ 



Ueardl = h(x) y 



(48) 



where h{x) is the frequency dependence (with x = 
hv/kTcMs), which simplifies to h{x) — —2 at the 
Rayleigh-Jeans part of the spectrum, y being called the 
comptonization parameter. The integral in the right-hand 
side of equation ( ^8| ) can be expressed as a function of 
as for the X-ray luminosity : 



AT 

TcMB 



= h{x)y 



h{x) 
h{x) 



kuT 

(Tt 
TOpC^ 



TieTdl 
(kTf^ 

K3/2 



dl 



(49) 
(50) 
(51) 



5.1. An analytic prediction for the y — T relation 



For current SZ experiments resolution (especially interfer- 
ometric ones), the beam smearing is still a critical issue 
as long as we want to determine a central value. Most 
authors use X-ray determined temperature as well as the 
X-ray surface brightness to correct for the beam and to 
obtain the central value of the SZ effect y{0). This might 
be a source of uncertainty but we will nevertheless express 
this quantity in our physical framework, and compare it 
to observations. 

The beam-corrected quantity can be obtained in the 
following way : going back to equation (|l]) we write (as- 
suming spherical symmetry) 



2/0 



2/(0) 
ijt 



7-5/2 



dl 



rdr 



(It 



-Rv 



dx 



Using eq. (|l]), ( |B.10| ) and (|1|) we have thus 



2/0 = 2/0 T'^ 
with 



2/0 



3M, 



-1 -3/2 



-,1/3 



4^Ap0r5' 
where we have defined : 
dx 



K, 



-3/2 



Q 



SZ 



7o P'^{^) 

The value of Qsz is computed in Appendix 



(52) 



(53) 



(54) 



(55) 



(56) 
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E 




Te (keV) 

Fig. 6. Central SZ decrement - X-r ay temperature rela - 
tion. The data points are taken from Zhang fc Wu (206o| ) . 
The soUd Hne is the predicted yo using Qsz = 1-835 . The 
dot-dashed hne corresponds to the renormahsed correla- 
tion fitted by Cooray (1999). The predicted normalisation 
was computed using Kq = 102keVcm^ and Tq = 2keV, 
i.e., it is not free once Kq has been fitted to the observa- 
tions (Sec. 2.3.3| ). 



Finally, let us make a comment about the shape of the 
y — T relation. For T ^ Tq, we get : 



yo oc oc t3/2 

while for T < Tq : 
yo (X T3. 



(57) 



(58) 



Thus, as a matter of fact we expect a steepening of the 
y — T relation when going from clusters to groups. 

5.2. Comparing with y — T observations 



Taking observationnal data from Zhang & Wu (200C) (di- 



amonds), we can compare our predicted relation of equa- 
tion (54) to observed correlations. 

In figure ^we plot data and their associated error bars, 
our prediction for Qsz — 1.835 (solid line, computed in 
Appendix]^) and the best fit of Cooray (1999, dot-dashed 
line)|^. The fact that only rich cluster data are yet avail- 
able (which limits the leverage on the slope) and the large 
dispersion on the observed y-T correlation translates into 
a large uncertainty on the best-fit slope : 2.35 ± 0.85. 

Our prediction for the beam-corrected value is steeper 
than the Cooray's best-fit relation but still consistent with 
it at the 1 a level. It seems to be in better agreement with 
Zhang's data set, which extends to slightly higher and 
lower temperatures than Cooray's. Note that we will fur- 
ther show that equation ( ^4|) is also in very good agree- 



Note that in jCooray (1999 )'s paper, the quoted normah- 
sation of the best-fit y — T relation is wrong by a factor of 
10. 



ment with semi-analytic models for both groups and clus- 
ters. 

The overall agreement with both data sets is less im- 
pressive than for the previously derived Lx — T relation, 
but is still very reasonnable, given the observational dis- 
persion. 

Since X-ray and SZ observations are completely inde- 
pendent, the agreement between our analytic model with 
the latter is another independent confirmation of its va- 
lidity. A robust confirmation would come from SZ ob- 
servations of groups of galaxies (unfortunately not yet 
feasible) and from the change in the correlation slopes. 
Nevertheless, we can compare our simple model to more 
elaborate semi-analytic models and their predictions. 

Menci & Cavaliere (2001) have presented semi-analytic 
predictions for the observed SZ effect in groups and clus- 
ters. Their semi-analytic model takes into account the pre- 
heating due to star formation, in order to bend the Lx — T 
relation in the groups mass range. Our predictions for the 
y{0 ) — T relation (namely equation ( |57| ) for clusters and 
( p8|) for groups) are in complete accord with their pre- 
heated models. They found a self-similar relation for high 
mass systems (y cx T^l"^) which is exactly our prediction. 
The slope in the groups mass range depends on the feed- 
back model they chose, but our prediction (j/(0) cx T^) is 
also similar to their computed moderate feedback model 
(which looks to be a better choice, since their strong feed- 
back model is likely to give only an upper limit on the 
feedback effect). We are thus able to reproduce their fig- 
ures 3 and 4 with our simple analytic model. Even if this 
is not a confirmation of the model (since we use an ana- 
lytic version of their model of shocks and our hypothesis 
on the preheating are very similar), we find it extremely 
satisfying to reproduce analytically, and in a physically 
straightforward and motivated way, the results of a more 
complicated and non-analytic model. This gives us more 
confidence in the fact, that we have captured in this sim- 
ple physical scheme the essential ingredients of the clusters 
and groups formation. 



6. Linking SZ and X-ray measurements 

To link SZ and X-ray measurements, it is natural to deter- 
mine a correlation between y and Lx ■ It is expected to be 
much tighter than the y-T relation (see Cooray, 1999), but 
is more intricate to obtain analytically. In this section, we 
will give two analytic expressions for this correlation in the 
rich clusters's and in the groups' range mass respectively. 
Indeed using equations ( ^^ and (|37|), we can write : 



vl 



(59) 



To proceed further we need now to "invert" equation (54) 
but this can not be easily done without a new hypotesis. 
Since the SZ observations span only the range of clusters 
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- SZ decrement (niK) 



Fig. 7. Central SZ decrement - X-ray luminosity relation. 
The solid line represents the analytic predicted correla- 
tion, using the predicted values of Qsz and Q x found 
in Appendix |c[ The dashed line corresponds to Cooray 
(1999|)ls best-fit correlation. The data are taken from 



Zhang fc Wu (2000) 



yet, as long as we want to compare to observations it is 
satisfying to assume that T ^ To. We can then write : 



r = T, 



-3/2 



2/3 



and inserting this in equation (p9) gives 



Lx = 



J a ^3/2 / 2/0 



_ ^0 .Tn3/2 



yo 



4/3 



^ with a =1/2. 



(60) 

(61) 
(62) 



The observed correlation between y and Lx, together 
with the analytic prediction (equation ( |62|) , solid line) is 
shown in figure The data are taken from [Zhang fc Wil 
(2 000|)]and the d ashed line is the best-fit relation found 
by Cooray (1999|) . The normalisation has been computed 
using Ko = 102keVcm2, Tq = 2keV, Qx = 0.115 and 
Qsz = 1.835. The agreement is very good, even if it must 
be recalled that the shown correlation has been computed 
for the case T ^ Tq. Thus the real correlation will be 
slightly steeper (thus more in agreement with Cooray's 
correlation), but by a very small amount, well within the 
observational dispersion. 

Finally, one can compute this correlation in the groups 
mass range (T <^ Tq), yielding a steeper correlation : 



Lx=L^xi'- 



3/2 



with a = 0. 



(63) 



7. Discussion 



We have shown that a simple modeling of the main process 
in charge of the formation of clusters (namely the shocks 



which thermalise the gas inside the virial radius) allows 
to derive a general K — T relation that reproduces very 
well the observed one. Moreover, we were able to deduce 
most of the observed scaling laws with the adequate nor- 
malisation, encompassing both groups and clusters. Thus, 
the physical meaning of e.g. the Lx — T and y ~ T rela- 
tions, and in particular of the break in their self-similarity 
when going from clusters to groups (only evidenced in X- 
rays for the moment) can be understood as a relic of the 
formation process of these systems : a competition be- 
tween an entropy floor prior to the shock and the shock 
process itself. This suggests naturally that entropy consti- 
tutes the best "observable" in clusters, a s already empha- 
sized by the modified spherical model of Tozzi fc Norman 



(2001 ) . What we add here is that we show that entropy (in 



fact, the adiabat) is also suitable to seek analytic expres- 
sions for the different correlations found observationally in 
groups and clusters. Reversing this argument, these scal- 
ing laws should be an appropriate probe of the entropy 
content of clusters. 



7.1. The Lx 
content 



T relation as a probe of entropy 



While the entropy spatial variations of the baryonic gas in 
clusters is a valuable piece of information on their present 
thermo dynamic state (in partic ular for merging clusters, 
see e.g., Markevitchet al., 1999 ), the total entropy content 
K (hereafter TEC), defined as 



K = I K{r)dV 
Jv 



(64) 



is related to the integrated thermodynamic history of 
the formation of these systems. Therefore, any non- 
gravitational entropy injection at some point in this his- 
tory will enhance the entropy content over the shock- 
generated one, while the cooling should decrease it in the 
central parts. 

The analytic model we have presented in the last sec- 
tions allows us to shed some light in the TEC and its link 
to the Lx—T relation. On the one hand, given our assump- 
tion of isothermality, eq. ( |3l| ) indicates that Lx/T^^" 
is related to the TEC in a given system. On the other 
hand, the self-similar models of cluster formation ( Kaiser J 



19861) predict that Lx ocT^, which gives directly (using 
a = 1/2) 



Lx 

y7/2 



oc 



rp-3/2 



(65) 



Any departure from this last relation witnesses entropy in- 
jection or loss. Consequently, at a given temperature, the 
gap between the curve representing eq. ( |65| ) and the locus 
of observed groups and clusters will give a indirect mea- 
sure of entropy gain at a given mass (if the central cooling 
parts are adequately excised). Unlike the observations of 
the "entropy fioor" , wh ich were obtained at a single scaled 
radius near the center ( Ponman et al., 1999 ; Lloyd-Davies 



|et al., 2000| ) and so are representative of a particular point 
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in the entropic history (i.e., at a time near the last major 
merger) , this measure includes most of the volume of the 
systems and thus most of its past thermodynamic history. 
However, this does not directly measure K, but the inte- 
gral of a power-law of the TEC. We will first show that 
in a semi-analytic scheme (with less assumptions than the 
present analytic model as well as with a realistic entropy 
evolution of the intergalactic medium) the interpretation 
of the Lx ~ T relation is still valid and then discuss how 
X-ray and SZ observations can constraint the TEC. 

7.2. Relaxing some hypotheses 

Assuming local thermodynamic equilibrium, isothermal- 
ity, and the self-similarity of the adiabat profile, we illus- 
trated how the observed Lx — T relation is a probe of the 
TEC. Whereas the LTE is quite robust (and actually very 
difficult in practice to alleviate), neither the isothermality 
nor the self-similarity are fully exact from both an obser- 
vational and a theoretical point of view. In particular, this 
last assumption is a very strong one and will be discussed 
in depth in section It is thus worth relaxing these 
hypotheses and testing the validity of our work by consid- 
ering how well the quantity K^^dV probes Ljf/T^/^, 
i. e. testing eq. ( |3l| ) when we consider a full non- isothermal 
and non-self-similar profile. 

To this purpose, we will use a semi-analytic model 
(hereafter SAM) developped independently by one of us 
(Dos Santos 2001, in preparation). This particular model 
uses the conduction-structured temperature profile (which 
was shown by Dos Santos (2001) to describe very well the 
temperature and surface brightness profile of clusters of 
galaxies) , together with an NFW profile for the dark mat- 
ter density profile. A shock model at the virial radius, 
together with entropic constraints at the center allow to 
predict the temperature profile without the hypothesis of 
isothermality or the unphysical polytropic link between 
temperature and density. The evolution of the central en- 
tropy is governed by the entropy evolution of the IGM 



1 .000 r 



in the universe obtained from the model of Valageas & 
Silk (1999D in two different cases : in the first one, the 



reheating is provided by SN explosions only, while in the 
second one, AGN and quasar activity provide the entropy 
injection. Both cases were validated against a number of 
observations, including the Lx — T relation, the change 
in surface brightness profiles from clusters to groups, the 
baryon fraction in these systems and the entropy floor. 

Figure ^ shows the comparison between observations 
and the SAM. The sofid (SN case) and dashed lines (QSO 
case) show the quantity K~^dV computed directly 
using the model. Since the specific entropy is known up 
to an additive constant, both lines were normalised to 
match the high temperature clusters. The dot-dashed line 
shows the self-similar prediction for Lx /T^^^ (eq. ^), also 
renormalised to match rich clusters. The preheated models 
high temperature slopes match naturally the self-similar 
prediction. They also match very well the cluster data 



0.010 - 



0.001 




T„ (keV) 

Fig. 8. Lx/T''/'^ (observational points) compared to the 
direct computation of jyK-'^(Y)dV in the SAM. The 
data are the same as in the figure |[ The solid (respec- 
tively, dashed) line is the result of the SAM in the case of 
SN (respectively. QSO) reheating. The dot-dashed fine is 
the self-similar prediction {Lx/T'^^'^ oc T~^/^). Thus, re- 
laxing the assumptions of isothermality and self-similarity 
of our analytic model does not change the interpretation 
of as a probe of JyK-^{v)dV. This explains 

that our simple model is able to recover the slope and the 
normalisation of the Lx — T relation. 



and their trend in temperature. It is in particular remark- 
able that the three clusters with the smallest temperatures 
(T ~ 2 keV) depart notably from the self-similar predic- 
tion, and lie exactly on top of the semi-analytic prediction. 
The lower envelope of groups is well followed by the SAM, 
while som e po ints are over this prediction. We think that, 
as in Sec. 



3.2 



this is due partly to the fact that the cen- 
tral cooling regions of groups were not removed, unlike the 
clusters, and partly to an intrinsic scatter. 

The figure shows why our simple model works well in 
recovering slope and normalisation of the Lx — T relation 



(and of the other relations we studied) : the equality (31 ) is 
also verified in a more general semi-analytic model where 
the adiabat profile is fully known and the ICM is not 
isothermal. At a given temperature, the difference in loca- 
tion between observations and the self-similar curve will 
give a measure of the non-gravitational entropy injection 
in systems with this mean temperature. This measure does 
not depend at all on hydrostatic equilibrium or spherical 
.symmetry. Nor does it require either that the preheating 
was external. In fact, eq. ( |3^ ) is completely independent 
of any modeling of the entropy injection. It only requires 
isothermality and local thermodynamic equilibrium on the 
scales probed by the observations. It thus can be used as a 
powerful, model-independent, constraint on all the inter- 
nal or external preheating models, as well as on differential 
galaxy forma tion efficiency between groups and clusters 
( Bryan, 2000| ). Unfortunately, the TEC is not measured 
directly, but the integral of the right-hand side of eq. (31) 
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will obviously decrease as the entropy injection amount is 
increased. Thus, if the adiabat profile is known in a given 
model, it will be easy to relate this integral to the TEC, 
even if a certain amount of degeneracy will obviously be 
present (this measure being integrated in space and time, 
different entropy injection histories can lead to the same 
final amount of entropy at z = 0). 

In fact, this constraint will be more or less the same 
as the constraint given by the Lx — T relation (already 
used to constrain the models), even if the interpretation 
in terms of entropy is physically more satisfying. However, 
the same remarks can be made with the surface integrated 
SZ decrement. Indeed, the SZ y-Compton parameter inte- 
grated over the whole surface of a cluster can be written : 



y 



(JT 



dV 




(66) 



which shows that the quantity y/T^/"^ is also a probe of the 
entropy content, but with a different power of the adiabat 
profile inside the integral. Thus, once SZ observations for 
groups are available, even non-resolved, the combination 
of X-rays and SZ at a given temperature will be a power- 
ful constraint on cluster formation models. Note that the 
relation Mjqm — T gives an analogous relation as eq. (|66|), 
and can thus be used now as an entropy probe. 

7.3. The adiabat profiles and the self-similarity 

To derive the analytic scaling relations presented in this 
paper, we have used a strong assumption : the self- 
similarity of rescaled adiabat profiles (see eq. ^2|) i.e., the 
fact that the temperature enters only the normalisation of 
the adiabat profile. The validity of this assumption is ques- 
tionable, both on observa tional and theoretical grounds. 

Ponman ct al. (1999) show that density profiles are 
shallower in groups than in clusters (even if this result is 



still uijcertain, see Roussei et ai., 2UU(J , tor an alternative 
view), while preheating models predict naturally that en- 
tropy gradients in lower temperature systems are smaller 
than in cl usters (due to adiaba tic infall during their life- 



T (keV) 



Fig. 9. Upper panel : Average scaled adiabat profiles for 
systems grouped by mean temperature (solid : 6-14 keV; 
dashed : 3.7-6 keV; dot-dashed : 1.3-3.7 keV; dotted : 
0.5-1.3 keV). Each individual profile was o btained from 
the best-fit gas and density profiles given by Lloyd-Davies 
et al. (200C| ) and scaled by [1 + T/Tq] /{1 + z).No trend is 



time, see Balogh et al., 1999; Tozzi fc Norman, 2001) 



However, these models predict also large temperature gra- 



dients in groups which are not observed (Tozzi & Norman, 



2001). Moreover, using ID hydrodynamic models with in- 
ternal or external preheating, Brighenti & Mathews (2001) 
obtain linearly rising entropy profiles in groups (outside 
an isentropic core for external preheating models) as well 
as in clusters. They thus produce naturally self-similar 
adiabat profiles in groups and clusters. The reasons of 
such discrepancies are unclear, but are certainly linked 
with widely differing simulation methods. Note however 



that Knight & Ponman (1997), using similar ID hydrody- 



namic HiimilaliouH, have found shallower slopes in groups 
compared Lu clusLers, wiUiouL Luking iuLu ucwuuL way pru- 
heating (whose effect is expected to widen the slope dif- 
ference) . 

Spatially-resolved spectroscopic observations of groups 
and clusters provide directly adiabat profiles in groups 



found with the temperature (unlike with the scaling used 
by Lloyd-Davies et al., 200C) and all the profiles are ap- 
proximately self-similar. The vertical dashed line shows 
the mean cooling-flow radius. Bottom panel : Scaled cen- 
tral entropy value as a function of the mean tempera- 
ture of each system. Again, no excess in groups is found, 
which validates our assumption of self-similar adiabat pro- 
files (at least in the core) and our approach. The dashed 
line shows the entropy fioor value infered from the data 
{Kq = 102keVcm2, corresponding to Ki = 120keVcm2). 



and clusters and can then help to settle this debate 
Unfortunately, the low surface-brightness of groups allows 
this type of study only to small radii^. Indeed, Lloyd 



Davies et al. (2000 ) obtained the adiabat profile for 12 
systems with T < 4keV up to r = 0.25 x r^. Even if this 



^ Note that this situation will be much improved with groups 
observations with XMM-Newton, due to its enhanced sensitiv- 
ity. 
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maximum extent radius is small, the study of the adia- 
bat profiles behaviour with temperature inside it is still 
worth, since the X-ray luminosity comes mainly from this 
central part of the cluster (due to the scaling of Lx)- 
These authors have rescaled their observed entropy pro- 
files by a factor T^^{1 + z), removing the effect of system 
mass and of the evolution of the mean density of the uni- 
verse. They found that the scaled adiabat profiles do not 
coincide, the less massive systems having higher scaled en- 
tropy profiles (see their figure 4). However, this particular 
scaling does not take into account the physical processes 
that change the adiabat normalisation with mass. What 
we have shown in this paper is that the adiabat profiles of 
groups and clusters should be rescaled {i.e., divided) by 
the quantity : 



1 



T 



(67) 



where the (1 + z)^^ factor comes from the redshift de- 
pendence of the temperature (T cx (1-1- z)) and the red- 
shift dependence of the central density (proportional to 
the mean density of the universe i.e., oc (1 -I- z)'^)|^. 

Using the best-fit /5— model for the gas density and 
the best-fit linear ramp for the temperature profile, we 
reproduced their adiabat profiles, and rescaled them indi- 
vidually by the quantity in eq. (|67|). Then we grouped the 
systems by t emperature and computed the mean scaled 
profile, as in Lloyd-Davies et al. (200C ). The scaled pro- 
files are displayed in t he upper panel of figure 7.3 (we used 
the same line styles as Lloyd-Davies et al., 200C ). No trend 
with temperature is visible now and the slight dispersion 
can be attributed to the assumption that the observation 
time equals the formation time of a system. Moreover, 
the four profiles are very similar between r ~ 0.05 x 
(the mean cooling-flow radius, inside which entropy loss 
is achieved by cooling) and the outer radius. This means 
that our assumption of scaled self-similar adiabat profiles 
is indeed right, at least in the core of the systems. To con- 
firm this, one can do the same scaling with the central 
adiabat values (at r = 0.1 x r„), which is shown in the 



bottom panel of figure 7.3. Instead of finding an excess for 
r < 4 keV as in the figure 5 of [Lloyd-Davies et al. (20001 ), 
we find now that all the systems lie on the same central 
adiabat symbolised by the dashed line and giving directly 
the value of Kq in our model {i.e., Kq — 102keVcm^). 

This shows that, when rescaled conveniently, the cen- 
tral adiabat profiles of groups and clusters are indeed self- 
similar and thus that our assumption is valid, as well as 
our whole approach of taking into account not only hy- 
drodynamic shocks but also an entropy floor established 
before them. This explains why our slope and normalisa- 
tion work so well for the Lx—T and ya—T relations, which 
put more weight on the core of the systems. It also enlights 
the fact that our shape factor for the gas mass is less accu- 



rate, since most of the mass lies at large radius where the 
self-similarity assumption breaks down (see Appendix 

7.4. Comparison with other results 

We have compared most of our analytic results with pre- 
dictions from semi-analytic models and found good agree- 
ment (in particular for the trends of the gas fraction with 
total mass and the SZ effect with temperature in low tem- 
perature systems). Let us now compare the value of the 
entropy floor we need to reproduce the curvature of the 
Lx — T relation with values assumed in other theoretical 
models of energy injection. Balogh et al. (1999| ) 



assume 



a constant entropy injection value of ~ 350 keVcm^ and 



Tozzi & Norman (2001 )'s model needs a value in the range 
190 — 960 keVcm"^ in order to steepen Lx — T. These val- 
ues, as noted by Lloyd-Davies et al. (2000), are higher 



than the observed value, which are likely to be upper 
limits. On the other hand, our model needs a value of 
Kq = 102 keVcm^, well within the observational range of 
70- 140/1^0^ keVcm^. 

As we have indeed fitted this value to the observa- 
tions, it may be asked if we do not force in fact this agree- 
ment. The answer is clearly no, as, first, the agreement of 
the K — T relation with the observations does not garan- 
tee that the same value of Kq will provide the a good 
description of the other relations (both their shape, via 
the temperature at which the self-similarity is broken. To, 
and their normalisation) which range from radio to X- 
ray data. Second, the other models reproduce as well the 
K ~ T relation, together with the Lx — T one (at least 
[Tozzi &: Norman, 2001 , 's model is compared to these ob- 
servations) with the same high value of the entropy floor. 
A lower value compatible with observations would not fit 
this relation. 

The main difference between our present work and 
these models is the fact that their low temperature sys- 
tems are entirely isentropic i.e., no shocks occur at all 
in groups. This is the claimed approximation made by 



Balogh et al. (1999), who want to investigate a limit- 



LIoyd-Davies et al. (2000 ) do not take into account the 
temperature redshift dependence since they divide by T. 



ing model (which naturally fails also in the clusters' mass 
range). However, Tozzi &: Norman (200 1[ ) modelise shocks 
and adiabatic infall, and it is not clear why shocks do not 
raise the adiabat profiles of their groups. 

This clearly shows that shocks can not be completely 
supressed in groups unless an unacceptably high entropy 
floor is needed in order to break the self-similarity in the 
Lx — T relation. Our model takes shocks and the entropy 
floor into account by construction both in clusters and 
groups and is thus able to reproduce nicely all the X-ray 
correlations with the observed value of Kq. 



8. Conclusion 

Throughout this paper, we have shown that the adiabat 
constituted a "key observable" in the ICM, not only be- 
cause it is a record of the thermodynamic past history 
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of diffuse baryons in groups and clusters, but also be- 
cause it allowed us to derive for the first time unifying 



injection level will affect the properties of a given group), 
but the number of groups detected will be large 



analytt c expression for the slope and normalisation of ob- 
served correlations of groups and clusters in X-rays and 
SZ. We have thus derived analytic expressions and their 
normalisation for the Lx ~ T, Mjcm — T, fgas ^T, yQ — T 
and yo — Lx correlations (for the latter, two limiting ex- 
pressions have been provided in the case of groups and of 
clusters respectively) . The Appendix ^ summarizes these 
correlations and gives the corresponding numerical values 
of their normalisations. We will conclude by discussing 
some applications of the model. 

Clusters of galaxies, and in particular their baryonic 
content observed through its X-ray emission, have often 



Rome: 



et al. (2001) have estimated that more than 100 clusters 



been used to provide cosmological constraints (Perrenod 
19801 ; [flCaiser, 1986| ; pukbir et al., 1997| ; [Blanchard et al. 
2000, i,nd references therein). These studies have shown 
that the cluster population can provide tight constraints, 
but have also revealed two problems (apart from purely 
statistical problems due to large error bars of e.g., the 
temperature in high-redshift clusters) : first, the difficulty 
to model the luminosity in a cluster, mostly due to the fact 
that it depends on the core properties, whose formation is 
still uncertain. Second, the small available number of rich 
clusters, in particular at high redshift, which precludes yet 
a reliable use of these test^. 

We think that the present work, together with the new 
generation of X-ray observatories, can alleviate both of 
these problems. First, our modelisation of the cluster lu- 
minosity does rely on the physical processes being at work 
during the formation of a cluster, and not on an hypothet- 
ical density profile. Therefore, the luminosity does not de- 
pend on an arbitrary core radius and is thus more reliable 
than previous attempts (in particular, we don't have to 
assume anything on the core radius evolution like e.g., 
Romer et al., 2001 ). Second, we extend the analytic re- 
lations to the groups mass scale, taking into account the 
competition between an external entropy floor and hy- 
drodynamic shocks (thus recovering the curvature of the 
Lx — T relation) and make predictions for other proper- 
ties of groups. This is useful since the group population 
is thought, in a hierarchical structure formation cosmol- 
ogy, to be much larger than rich clusters. However, its 
redshift evolution has never been used to constrain cos- 
mology, partly because of the lack of reliable data and 
partly because no analytic modeling of the groups bary- 
onic properties was available (extrapolating the cluster 
baryonic properties was obviously wrong, as the recent 
observations show). The first problem will be solved, at 
least locally, by Chandra and especially XMM-Newton. 
Our model provides an attempt to solve the second prob- 
lem. Of course, groups will certainly be much more dis- 
persed than clusters (because the entropy injection will 
have more effect on them, and any spatial variation of this 



It is paradoxical but amusing to note that this same small 
number of high redshift cluste rs becomes a powerfu l test of the 



total density of the universe (Donahue et al., 1998) 



with T > 2 keV will be detected in a serendipitous sur- 
vey for z < 0.2 (irrespective of the cosmology), a number 
10 times highers than for cluster with T > 4keV and 50 
times highers than for T > 6 keV for a survey surface of 
- 800 deg^ The number of groups with T < 2 keV will ob- 
viously be superior, certainly not by a factor 10 because 
of the reduced size and steepening of the Lx — T rela- 
tion, but at least by a factor of a few. This number of 
well observed groups (together with pointed observations) 
will allow an unprecedent local calibration of the group 
ensemble properties. Using our model, it will be possible 
to use groups and clusters to constrain cosmology and the 
cluster LF, TF and Lx — T evolution to add further con- 
straints. Moreover, groups will provide useful constraints 
on the total amount of reheating the present universe un- 
derwent and its redshift evolution if high redshift groups 
are available (see Sec. 0). 
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Appendix A: Numerical values for the correlations 

In the following appendix, we provide numerical expres- 
sions for the normalisations of the different correlations 
the present model predicts. These will thus be straight- 
forward to use in another context. The temperature T is 
in keV, and the normalisations were computed using a lo- 
cal Hubble constant of Hq = 100 /i2/3 kms""'^ Mpc~^ with 
^2/3 = 2/3. We included the dependence on /i2/3 in each 
of the following expression^ 

We first describe one relation (external to the model 
and independent from it), the M — T relation, that we 
used to relate the virial radius to the temperature. Then, 
we write the different correlations presented in this work, 
namely Lx - T, Mjcm - T, fgas -T,yo-T and yo-Lx- 

A.l. The Mass-Temperature relation 

This relation is crucial, in that it allows us to relate the 
virial radius (where the shock is assumed to take place) 
with the temperature of the system, via the definition of 
M in terms of . The simulation results of ( [Evrard et al. J 
1996) are frequently used for mass-temperature scaling, 
but they seem to provide systematically a higher normal- 
isation than the observed one, whatever the method used 
to measure M (Nevalainen et al., 200C, and references 



therein). We used the observational results of Nevalainen 
[et al. (2000| ), derived from observed density and temper- 
ature profile, because their sample, although small, goes 
all the way from groups to clusters. These authors give 
the M — T relation at different scaled radii, from r2ooo 
to rsoo (where ta is the radius whose mean interior den- 



sity is A times the critical density). Since Evrard et al 
(1996D have shown that, inside a radius rsoo, the baryonic 



gas is in hydrostatic equilibrium to a very good approxi- 
mation and that hydrostatic masses measured within this 
radius should be reliable, we choose this radius to nor- 
mali se the M — T relation. M oreover, while at smaller ra- 
dius Nevalainen et al. (2000| ) find that the M — T relation 
is significantly steeper (99.99% confidence at riooo) than 
the self-similar prediction (M cx T'^/^), the slope they ex- 
trapolate to rsoo is consistent with 3/2 {x^ /f = 5.1/5, but 
note that the best fit has a slope of 1.84). We thus used 



* We thus have /150 — 4/3 and /iioo 
hso = iHo/50) and hioo = (Ho/100). 



2/3 where 
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throughout the paper the foUowing relation 

3/2 

M50O 



M.I J 



= 9.6 X IQi'^Mo/i^/l I — 



T 



3/2 



(A.l) 



which defines the values Af* and used in the text. Each 
time we speak about the virial radius in this paper, we 
refer to the radius rsoo related to M500 defined above, and 
we will thus take A — 500 throughout. 

Finally, let us state that introducing a steeper slope for 
the M — T relation is straightforward within the physical 
scheme we used in this work, but does not change much 
the predicted relations gathered in this appendix and their 
agreement with observations. Since we wanted to discuss 
the generic effect of shocks and preheating, the same ef- 
fects which are thought to steepen the M — T relation, we 
found natural to keep a self-similar scaling. 

A. 2. The predicted X-SZ-T correlations 

We now compute the different normalisations given in 
eq. @, @, @ and @. For this purpose, we 

will use the values of , and A from the last section. 
The critical density is defined as usual as : 



Ml 

SttG 



(A.2) 



and, for the luminosity, we will need the normalisation of 



the cooling function, taken from Eke et al. (1998) 
Ae(r) = 



= 1.2 X 10-24 J^l/2 



erg cm^ s ^ 



(A.3) 



where T is in keV. 

Finally, we take for To the value computed from Kq 
in the Sec. 2.3.2 i.e., Tq = 2keV and the values of the 
shape factors Qx, Qsz and Qm computed in Appendix 
Therefore, the only adjusted quantity used to produce these 
relations is Kq, which was fitted to the data of Hclsdon| 
fc Ponman (2000D in Sec. ^.3.3| to find = 102 kcVcm^. 
Once this quantity was fixed by comparison with the cen- 
tral entropy observations, no renormalisation of the fol- 
lowing relations is allowed. 

All the numbers are computed here at z = 0, but it is 
easy, within a given cosmology, to extend the formulae in 
the text to higher redshift (note that the value of Ko{z) 
must be assumed as well). The analytic correlations found 
are then 



• Lx — T relation 



Lx = 5.643 X 10"*^ 



• MicM T relation 



-to 



1 -3 



/i2/3ergs 1 (A.4) 



MicM = 2.405 X 10'^ T-' 



T 
Jo 



-1 -3/2 



/i-/3 Me (A.5) 



• fgas T relation 



- 0.0792 r3/2 



1 



T 



-3/2 



'2/3 



• yo — T relation : 

yo = 3.506 X lO^^^T^ 



1 



T 



-3/2 



-2/3 
''2/3 



(A.6) 



(A.7) 



which can be translated in central temperature decre- 
ment as (using Tcmb = 2.73 K and the RJ approxi- 
mation) 

, -| -3/2 

AT(0) = -1.914x10-2^3 ■ 



-to 



hyl^mK.{k.\ 



• AT(0) - Lx relation : 

★ for T ^ To (rich clusters temperature range) : 



L 



X 



2.206x10 



45 



-Ar(o)V/^ _io/9 



erg s 



ImK ; 2/3 
★ for T <C To (small groups temperature range) : 



i(A.9) 



Lx = 2.133x10 



45 



-Ar(o) \ 

ImK ) 



3/2 



/i2/3 erg s-^(A.lO) 



Appendix B: Infall velocity and mean temperature 

In this Appendix, we will find an expression for the infall 
velocity before the shock (noted vi) as a function of the 
mean temperature of the accreting system. To this pur- 
pose, the following reasonning in two steps is appropriate : 
we first show that v\ oc 0^ provides a good approximation 
and then derive a usefuU expression for 
We will then compare our result to hydrodynamic simula- 
tions, which will validate the use of the analytic expression 
found in the groups mass range. 



B.l. Reexpressing tlie in falling velocity 

Let's assume that the gas inflow is stationary. In a spher- 
ical collapse model, the gas is assumed to be at rest at 
a radius rta (the so-called "turn-around radius") before 
falling into the cluster potential. Since the gas will not be 
subjected to many processes changing its internal energy 
during this inflow, we can assume that the flow is isen- 
tropic. Thus, applying the Bernouilli equation between rta 



and ry (just befo re the shock) we obtain (see e.g. Landat 
Lifshitz, I959I) : 



.2 



+ UJ1+ (f)u = UJta + 



(B.l) 



where uj denotes the gas specific enthalpy (a; = CpT with Cp 
the specific heat at constant pressure). Denoting pta and 
pi the densities at rta and r.^ and , the following expression 
for vi can easily be found : 

2/3N 



-2(j)v 



1 



Ti 



c„— I 1 - ( ^ 

Pi 



Pta 



(B.2) 
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In the spherical collapse model, the density at r„ is 
times larger than at rta- We can thus write : 



'2(l)y 



1 



4 6„ 



(B.3) 



Obviously, in large mass systems, the last ratio will be 
negligible i.e., the thermal energy of the gas will be negli- 
gible when compared to its potential energy (the so-called 
"cold inflow" hypothesis). As gas is preheated {i.e., as Ti 
rises) or as the mass of the system is lowered this assump- 
tion is questioned, since the thermal content of the gas 
will be of the same order of the kinetic and potential en- 
ergy. Nevertheless, we will assume in the following that 
this limiting behaviour only occurs at mass scales smaller 
than groups of galaxies. We will show at the end of the ap- 
pendix that hydrodynamic numerical simulations of struc- 
ture formation in different cosmological models validate 
this last assumption. 
We have thus : 



(B.4) 



To take into account some uncertainties and shortcomings 
of this crude treatment (as the fact that the infall is usu- 
ally along filaments and so not spherical), we will write 
the slightly more general formula : 



where rj is considered as a constant, rj could be calibrated 
against numerical 



Miniati et al. 



simulations (see e.g., 
2000|)~^r chosen so that the flow at inflnite is the Hubble 
flow, as in Cavaliere et al. (1998| ). When used to compute 
the value of To in Sec. 2.3.2| , 77 is determined by the former 
method. 

Note that the established fact that the square of 
the infall velocity is proportionnal to the gravitational 
potential, both taken at the virial radius, is actually not 
a surprise. Indeed, before the shock, the gas is thought 
to follow the dark matter evolution, its velocity being 
equal to the dark matter on e up to the virial r adius (see 
for example the figure 15 of Frenk et al., 1999 ). Since it 
has been shown that the dark matter infall velocity in a 
spherical model scales as M^^^, while the mea n potential 
(j) cx M/R sc ales as M^/^ (|Ryu & Kang (1997|) , based on 
the work of Bcrtschinger (1985 )), it is thus natural that 



cx Af cx (/) 



B.2. Linking (j)^ with T 

We now look for a useful expression for without as- 
suming any particular analytic expression for (/'(r). First, 
since the Poisson equation states : 

V^,/) = 47rG'p, (B.6) 
we can write any gravitational potential in the form : 
<l>{x) = AnGp,rl4>{x), (B.7) 



where pc is a characteristic density (the background cos- 
mological density for example), is the virial radius, G 
the gravitational constant and (pi^) is a dimensionless po- 
tential with X = r/vy. Therefore we have in particular. 



='iTTGpcrl^{l) 



(B.8 



By using this as a general relation, we assumed implicitely 
that (^(1) does not depend on the mass of the system. 
However, numerical simulations indicate that the general 
expression for the gravitational potential is indeed uni- 
versal {i.e. it takes the same form in different cosmologi- 
cal models), but its normalisation should depend slightly 
on the mass via the so-called concentration parameter c 



(Navarro et al., 1997, hereafter NFW). Indeed, using the 
so-called NFW ana lytic expression, (p^ can be expressed 
as (see for example Lokas fc Mamon, 2001 ) ^ : 



= "^T^Gpc 



lie) 



(B.9) 



Since in these models, c is a slowly varying function of 
the mass (c(M) = (M/Mi5)-*' 086 for a ACDM universe 
at z = 0, M15 being the mass in units of 10^^ Mq) it is 
easily shown that g{c) cx T" oi2^ ^^ig^^ jg gg^y^ ^^^^ ^^^^^^ 
be considered as constant to a very good approximation. 

As a consequence^ ignoring the mass dependency of 

(^(1) in equati on ( B.g) is reasonnable and we will thus 
use equation (B. 



as a definition of the virial radius 



(B.5) gravitational potential. 



An even more practical form for equation ( |B.q ) does 
involve the (supposed isothermal) ICM temperature T. 
It is achievable through the use of the virial theorem. It 
indeed provides us with general relations between ry — M 
and Af — T in the following form : 



3Af 



in Apc 



1/3 



M = 



3/2 



(B.IO) 
(B.ll) 



where A is the virial radius overdensity (taken through- 
out this paper to be 500, see Appendix |^), pc is the 
critical density and and are constants defined in 
Appendix As a conclusion, we can now write equa- 
tion (p!8|) as : 



(j)y = ATrGpc<p{l) 
T 



ilT Apc 



2/3 



(B.12) 
(B.13) 



Eventually, by combining equations ( B.13| ) and (B.5) 
we reach the following relation describing the scaling of v1 
with T : 



= -27701 



(B.14) 



' Note that t he function q(c) defined here differs from the 
one adopted by Lokas & Mamon (2001) 
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Equation (B.14) is particularly interesting since it can 
be directly compared to the the results of hydrodynamic 
numerical simulations. Miniati et al. (200C ) have examined 
the properties of shock waves around clusters and groups 
of galaxies in two hydro+N-body realisations of two dif- 
ferent cosmological models. Having selected all structures 
with X-ray luminosity greater than 10"*^ erg/s, they found 
that the relation between the mean temperature and the 
infall velocity is given by an equation of the form (we only 
quote their result for the ACDM cosmology, analogous ex- 
pressions being found in the SCDM cosmology) : 



0.52 



(B.15) 



with 



VsO 



1.9 X lO^kms"^ and 



Ts = 7.8 X 10' K 

Their result spans a range from 10^ to several 10'' K, i.e. 
from poor groups to clusters, and is analogous to equa- 
tion ( B.14| )F^. This, together with the fact that secondary 
infall models predict the same type of relation as the one 
we deduced between vi and validates the assumption 
of cold inflow that was made in deriving equation 
even for the group mass range. 



Appendix C: The shape factors Qx, Qsz and Qm- 

In this appendix, we evaluate the shape factors Qxj Qsz 
and Qm, which enter the normalisation of the Lx — T, 
y — T and Mjcm — T relations. These quantities could in 
fact be considered as renormalisation constants obtained 
by comparison with observations. We will instead show 
here that they can be computed using a motivated model 
for the adiabat profile. The model we present thus predicts 
not only the shape but also the normalisation of the above 
correlations. 

C.l. Normalisation temperature 

The normalisations of the three relations cited above are 
defined in equations (|38|), (42) and (p^), where the tem- 
perature is in keV. The temperature at which the nor- 
malisation must be evaluated is obtained by equating the 
temperature shape functio n in each of these relations to 1. 
This is illustrated in figure CT , which shows the following 
functions (the last two ones being the same) : 



Lx(T) = T (solid line) 

[1 + T/Tof 



y{T) 



r3 



(dashed line) 



[l + T/Taf^ 
Micm(T) = ,3/2 (dashed Ime), 



[1 + T/To 



T (keV) 



Fig. C.l. Normalisation temperature of the Lx — T (solid 
line), y — T and Mjcm — T (dashed line) relations. The 
three relations will thus be normalised at T ~ 1.3 keV, in 
the groups mass range. 



and their intersection with the horizontal unit line. We 
can see that the normalisation temperature of the three 
relations occurs at T 1.3 keV. These temperatures are 
well into the groups temperature range, and we will thus 
compute the normalisation in this regime. First of all, we 
need an motivated analytic adiabat profile to compute the 
integrals defining the three normalisation constants. 

C.2. An analytic model for the adiabat Profile 

Groups have a low surface brightness in X-rays and are 
thus difficult to observe up to their virial radius. Obtaining 
a temperature profile at this radius is obviously even more 
challenging. While the new generation of X-ray satellites 
(Chandra and XMM-Newton in particular with its excep- 
tional sensitivity) is expected to produce constraints on 
low density regions, ROSAT (whose energy band was per- 
fectly adapted to the groups observation) has provided 
data mostly on the inner pa rts of these systems (see 



Helsdon fc Ponman, 200"c| ). Nevertheless, [David et al 



1996|) have published adiabat profiles for two cool clus 



ters (T ~ 1.5 keV) nearly reaching the virial radius (as 
well as for even cooler groups, see their figure 10), using a 
/3— model for the density and the best-fit power-law model 
for the temperature (despite large error bars). They found 
a isentropic core, followed by a very modest increase com- 
pared to clusters (where the virial value of the adiabat is 
at least 20 times higher than the core value). 

This can be easily understood, if one recalls that 
groups are thought to accrete most of their gas adia- 
batically, and should thus have constant adia bat profiles 



( Palogh et al., 1999| ; [Tozzi fc Norman, 20011) . However 



We will take the exponent of eq. ( ^3.15[ ) to be 1/2. 



this produces unacceptably large temperature gradients 
in small groups and the quasi-absence of the co oling- 
flow phenomenon , except in the largest clusters (Tozzi 
|fc Norman, 200 1" , see their figures 4c and 7b). It shows 
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that somewhere between the center and the virial radius, 
shocks must increase somewhat the adiabat, in order to 
lower the temperature gradient, giving rise to an adiabatic 
inner core (much larger than the core radius) followed by 
a modestly increasing adiabat profile]^. 

One can seek an analytic expression for the adia- 
bat profile by making use of the isothermal /3— model 
(Cavaliere & Fusco-Femiano, 1976). Suppose that the gas 
density profile is describe by this model i.e., 



n{x) = riQ 



l+(^ 

Xo 



21 "3/3/2 



(C.l) 



where the radii are rescaled in units of and xq is the core 
radius. Assuming isothermality, the dimensionless adiabat 
profile is then 



fix) 



Xq 



(C.2) 



Using the semi-analytic model already discussed in the 
text (Dos Santos 2001, in preparation, see Sec. [z^ ), we 
found that, in an external preheated model (by QSOs or 
SNs) , the adiabat profile of groups and clusters can be de- 
scribed approximately by the last equation with (3=1/2 
and a core radius varying with mass. Even if the systems 
modelled are not isothermal, this result is not surprising as 
the shocks are expected to leave a linearly rising adiabat 
profile (hence the value of (3, see Thomas & Couchman 
1992 |;[|ozzi &: Norman, 200l| ; [Brighenti fc Mathews, 200l| ) 



while the isentropic accretion of gas leaves a constant adi- 
abat profile. Fixing /? to 1/2, the core radius xq deter- 
mines the extend of the central isentropic core. We have 
found that Xq = 0.7 reproduces approximately the adiabat 
profile of small groups in our semi-analytic model, while 
X = 0.1 reproduces its behaviour in rich clusters. This is in 



line with the observational results of David ct al. (1996) 
showing that the profiles are much shallower in groups 
than i n clu sters. Thus we will use = 1/2 and xq = 0.7 
in eq. ( |C.2| ) to compute the values of the shape factors. 

C.3. The shape factor Qx 

The theoretical shape factor is defined by: 



Q 



X 



X 



^dx 



JO 
'0 

0.1150 



1 



f-V 

V0.7/ 



-3/2 



This value is very close to the best-fit value (0.12, consider- 
ing Qx as an adjustable parameter) and undistinguishable 
given the observational errors. Although the normalisation 

Note that HCG62, one of the best examples of relaxed 
group (T ~ 1 keV) , has a ste ep declining temperature out - 



side the inner core (see fig. 4 of Finoguenov fc Ponman, 1999|) , 
but has also an obvious cooling-fiow. 



has been computed in the groups mass ranges, it recovers 
perfectly the slope and normalisation of rich clusters. We 
will thus use Qx = 0.115. 

C.4. The shape factor Qsz 

The shape factor for the y — T relation is given by : 



Qsz 



/•OO 








1 + 


(^)l 


lo 





-3/4 



dx 



(C.3) 



1.835 



Here again the normalisation is in very good agreement 
with the observational data from Zhang fc Wu (2000 ), 
despite the large dispersion of the latter (see figure ^ . 



C.5. The shape factor Qm 

The shape factor for the Micm — T relation is given by : 

2-1 -3/4 



Qm = I X- 
/o 

0.1915 



1 



X N 2 

OJ 



The best-fit value to the data of [Mohr et al. (1999| , here- 
after MME) is Qm = 0.11 i.e., ~ 40% lower than the 
value we compute. However, apart from pure measure- 
ments uncertainties and known biases^, several system- 
atic effects can cause this disagreement : first, the major 
part of the 45 clusters they study are not detected up to 
the radius tsoq. Thus, they extrapolate the data taking 
their best-fit slope out to rsoo, which can be a source of 
errors. Second, as MME explain, computing r5oo requires 
a model of the potential or a rescaling to the value of a 
given cluster (they chose A1795 for the scaling of their 
virial relations). If one of the assumptions made is erro- 
neous, it is possible that their numbers quote a radius ta 
with A ^ 500. Indeed, if rsoo is underestimated by only 
20 %, we obtain : 



f-V 

V0.7/ 



-1-3/4 



0.1136, 



which is very close to the best-fit value. We use this last 
value for Qm in the Mjcm — T relation, as well as in the 
fgas - T relation. 

The luminosity of a cluster is dominated by the core 
gas (since it scales as the square of the density) and the 
SZ efi'ect computed here is a central value (taking into 
account the beam smearing effect), while the gas mass 
depends mainly on the outer parts of the profile, where 
ourassumption of self-similarity of the scaled ad iab at pro- 
files is more likely to break down (see section 7.3). This 



Note that testing their analysis method with hydrodynamic 
numerical simulations, MME find that they are likely to over- 
estimate the ICM mass by a factor of ~ 12 %, which increases 
the disagreement with our Qm value by the same amount 
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explains naturally why the predictions for Qx and Qsz 
are much more accurate than for Qm- Nevertheless, we 
consider the agreement between the theoretical values and 
the observed normalisations very satisfying. The fact that 
with the same model for the adiabat profile we can pre- 
dict the correct normalisations of the cluster correlations 
in two very different wavebands is obviously a sign of co- 
herence of the whole scheme. Since the shape factors are 
predicted, the only parameter of the model is the central 
entropy in groups Kq, which was fit to the observations of 
central adiabat value. 



